核心任务:
- 理解理论: 解释常微分方程解的“生命周期”(即能延伸多远——延拓定理)以及解对“起跑线”(初值)的敏感程度(连续性与可微性)。
- 应用实战: 利用定理5中的微分公式(公式3.17和3.18)解决两个具体的计算问题。
关键难点与隐含假设:
- 难点: 理解 ∂y0∂φ 这种符号的物理意义。它代表的是“当初值 y0 变化一点点时,整个函数曲线在 x 处的变化率”。
- 假设: 在计算题目时,必须先知道“原方程的解 φ(x) 是什么”,才能将其代入积分公式中进行计算。特别是第2题,隐含了一个特解 y=0。
直观理解:
想象你在一个大广场(区域 G)上按照一张藏宝图(方程 y′=f(x,y) 指示的方向)行走。
- 定理的核心含义: 只要你还在广场内,且指示方向是清晰的(f 连续且满足利普希茨条件),你就不会突然“卡住”或“消失”。你一定能一直走下去。
- 停下来的原因只有两个:
- 撞墙了: 你走到了广场的边缘(趋于边界)。
- 飞出去了: 你的速度变得无穷大,或者你走到了无穷远的地方(y 无界或 x 无界)。
例1的深度解读:
- 方程 y′=(y2−1)/2。
- 解 y=1+ex1−ex (过原点): 这个解被限制在 −1 到 1 之间,它永远不会变成无穷大,所以它可以一直向左右延伸到无穷远。这就是推论中的情况1。
- 解 y=1−ex1+ex (过 (ln2,−3)): 当 x→0 时,分母趋近于0,整个 y 趋向于无穷大。这意味着解“爆炸”了。这就是推论中的情况2(y 无界)。
直观理解:
想象射击。
- 连续性 (定理4): 你的手抖了一下(初值变动一点点),子弹的落点只会偏一点点,不会突然飞到背后去。
- 可微性 (定理5): 我们不仅想知道偏了,还想知道偏了多少。
- ∂y0∂φ:不仅是连续的,还是光滑的,我们可以精确计算出“初值改变 1 单位,解在 x 处改变多少单位”。
公式解读(至关重要):
我们要计算的是解 φ 对初值 x0,y0 的偏导数。这里有两个现成的公式(无需死记硬背,但要会用):
- 公式 (3.18) 对 y0 的偏导:∂y0∂φ=exp(∫x0x∂y∂f(t,φ(t))dt)
- 含义: 这是一个指数累积的过程。∂y∂f 就像是一个“放大系数”,如果它大于0,初值的微小误差会随着 x 增大而指数级放大。
- 公式 (3.17) 对 x0 的偏导:∂x0∂φ=−f(x0,y0)⋅∂y0∂φ
- 含义: 这个公式很有趣,它其实告诉我们,对 x0 的导数可以通过对 y0 的导数算出来。它多乘了一个 −f(x0,y0)(即初始处的斜率)。
我将一步步为你讲解这两道题目。
题目: 假设 P(x),Q(x) 连续,y=φ(x,x0,y0) 是方程 dxdy=P(x)y+Q(x) 的解。求 ∂x0∂φ,∂y0∂φ 及 ∂x∂φ。
详细解析:
Step 1: 识别 f(x,y) 及其偏导数
根据方程 dxdy=P(x)y+Q(x),我们可以确定右端函数为:
f(x,y)=P(x)y+Q(x)
我们需要计算 f 对 y 的偏导数 ∂y∂f,这是使用定理5公式的关键:
∂y∂f=∂y∂[P(x)y+Q(x)]=P(x)
Step 2: 计算 ∂y0∂φ (利用公式 3.18)
根据公式:
∂y0∂φ=exp(∫x0x∂y∂fdt)
将 ∂y∂f=P(x) 代入(注意积分变量用 t 以示区别):
∂y0∂φ=exp(∫x0xP(t)dt)(解释:这就是如果你的初始 y 值改变一点点,解在 x 处变化的大小。因为它是一阶线性方程,这个放大因子只与 P(x) 有关。)
Step 3: 计算 ∂x0∂φ (利用公式 3.17)
根据公式:
∂x0∂φ=−f(x0,y0)⋅exp(∫x0x∂y∂fdt)
我们需要先写出 f(x0,y0)。将 x0,y0 代入 f(x,y) 的表达式:
f(x0,y0)=P(x0)y0+Q(x0)
结合 Step 2 的结果,得到:
∂x0∂φ=−[P(x0)y0+Q(x0)]⋅exp(∫x0xP(t)dt)
Step 4: 计算 ∂x∂φ
这个最简单,但容易想复杂。φ 是原方程的解。
解的定义就是:把它代入方程,方程成立。
即 ∂x∂φ 就是 y′,也就是方程的右端:
∂x∂φ=f(x,φ)=P(x)φ(x,x0,y0)+Q(x)
题目: 给定方程 dxdy=sinxy,求 ∂x0∂y 和 ∂y0∂y 在 x0=1,y0=0 时的表达式。
详细解析:
关键提示: 对于非线性方程,积分里的 ∂y∂f 通常含有 y。要计算积分,我们必须先知道具体的解 y=φ(x) 是什么。
Step 1: 寻找通过点 (1,0) 的特解
我们观察方程 dxdy=sinxy 和初值 y(1)=0。
让我们猜测一个最简单的解:恒为 0 的解,即 y≡0。
- 左边 y′=(0)′=0。
- 右边 sin(0/x)=sin(0)=0。
- 左边 = 右边,且满足初值 y(1)=0。
因此,此初值问题的解就是 y=φ(x)=0。
(这一步非常关键!如果不知道解是0,后面的积分就没法算。)
Step 2: 计算 ∂y∂f 并代入特解
由 f(x,y)=sin(xy),根据链式法则求导:
∂y∂f=cos(xy)⋅∂y∂(xy)=x1cos(xy)
现在,我们将具体的解 y=0 代入上式:
∂y∂fy=0=x1cos(0)=x1⋅1=x1
Step 3: 计算 ∂y0∂y
利用公式 (3.18),注意下限 x0=1:
∂y0∂y=exp(∫1x∂y∂fdt)
将 Step 2 算出的 t1 代入积分:
∫1xt1dt=ln∣t∣1x=ln∣x∣−ln1=lnx(考虑 x>0)
所以:
∂y0∂y=exp(lnx)=x结论: ∂y0∂y=x
Step 4: 计算 ∂x0∂y
利用公式 (3.17):
∂x0∂y=−f(x0,y0)⋅∂y0∂y
我们需要计算 f(x0,y0),也就是 f(1,0):
f(1,0)=sin(10)=0
由于前面的系数是 0,无论后面是多少,结果都是 0。
∂x0∂y=−0⋅x=0结论: ∂x0∂y=0
- 理论核心: 只要记住,解对初值的偏导数 ∂y0∂φ 本质上是一个解的“变分”,它服从一个线性的微分关系,所以结果总是一个指数积分形式 (exp(∫…))。
- 做题技巧:
- 如果是线性方程(如题1),直接代入 P(x) 即可。
- 如果是非线性方程(如题2),必须先找出原方程的特解,代入 ∂y∂f 把它变成只含 x 的函数,才能进行积分。
这些题目覆盖了线性微分方程的四大支柱:
- 叠加原理(解的性质)
- 常数变易法(解的计算,核心中的核心)
- 初值问题(归一化)
- 理论证明(朗斯基行列式与解的空间结构)
核心思想:
微分方程其实是一个“机器”(算子)。你把函数放进去,它就吐出一个结果。
这道题是说:如果把 x1 放进去得到 f1,把 x2 放进去得到 f2,那么把 x1+x2 一起放进去,得到的就是 f1+f2。
证明步骤:
第一步:简化符号(为了不晕车)
为了不每次都写那一长串 dtndnx+⋯,我们定义一个记号 L[x]。
令 L 代表方程左边的运算:
L[x]=dtndnx+a1(t)dtn−1dn−1x+⋯+an(t)x
第二步:利用导数的“线性”性质
我们要知道一个常识:和的导数 = 导数的和。
即 (x1+x2)′=x1′+x2′,(x1+x2)′′=x1′′+x2′′,以此类推。
这意味着运算 L 也是线性的:
L[x1+x2]=L[x1]+L[x2]
第三步:代入题目条件
题目告诉我们:
- x1 是第一个方程的解 ⟹L[x1]=f1(t)
- x2 是第二个方程的解 ⟹L[x2]=f2(t)
第四步:得出结论
我们要验证 x1+x2 是不是第三个方程的解,就是把 x1+x2 塞进 L 里看看:
L[x1+x2]=f1(t)L[x1]+f2(t)L[x2]=f1(t)+f2(t)得证。
利用常数变易法求解。
Step 1: 检查系数x′′ 前面系数是 1,合格。f(t)=cost。
Step 2: 列方程组(常数变易法模板)
我们需要求 c1′,c2′。
[
{c1′et+c2′e−t=0c1′(et)′+c2′(e−t)′=cost⟹{c1′et+c2′e−t=0—(1)c1′et−c2′e−t=cost—(2)
]
Step 3: 解方程组
(1) + (2) 得:2c1′et=cost⟹c1′=21e−tcost
(1) - (2) 得:2c2′e−t=−cost⟹c2′=−21etcost
Step 4: 积分(这里考微积分基本功)
这两个积分需要用到“分部积分法”或者“欧拉公式”。
(注:∫eatcostdt 是经典循环积分,公式为 a2+1eat(acost+sint))
- c1=21∫e−tcostdt=21[2e−t(−cost+sint)]=41e−t(sint−cost)
- c2=−21∫etcostdt=−21[2et(cost+sint)]=−41et(sint+cost)
Step 5: 写通解
通解 x=C1et+C2e−t+(c1x1+c2x2)。
特解部分:
c1et=41(sint−cost)c2e−t=−41(sint+cost)
特解 xˉ=41sint−41cost−41sint−41cost=−21cost
最终结果: x=C1et+C2e−t−21cost
Step 1: 检查系数x′′ 前系数为 1,合格。f(t)=tsin2t。
Step 2: 列方程组
[
{c1′cos2t+c2′sin2t=0c1′(−2sin2t)+c2′(2cos2t)=tsin2t
]
Step 3: 解方程组
这是一个二元一次方程组。
- 将第一式乘以 2sin2t,第二式乘以 cos2t,相加消去 c1′:
c2′(2sin22t+2cos22t)=(tsin2t)(cos2t)2c2′=tsin2tcos2t=21tsin4t⟹c2′=41tsin4t
- 求 c1′:由第一式 c1′=−c2′cos2tsin2t=−41tsin4tcos2tsin2t=−21tsin22t
(利用三角公式降幂:sin22t=21−cos4t)
⟹c1′=−41t(1−cos4t)=−41t+41tcos4t
Step 4: 积分(略去繁琐步骤,你需要自己练习分部积分)
这里积出来会比较长,但逻辑就是简单的积分。
Step 1: 检查系数(陷阱!!)x′′ 前面是 t2,必须两边同时除以 t2!
标准化方程:
x′′−t1x′+t21x=t6+34
所以,f(t)=t6+34。如果你用 6t+34t2 算,就全错了。
Step 2: 列方程组x1=t⟹x1′=1x2=tlnt⟹x2′=lnt+1
[
{c1′t+c2′tlnt=0c1′(1)+c2′(lnt+1)=t6+34
]
Step 3: 解方程组
由第一式:c1′=−c2′lnt。
代入第二式:
−c2′lnt+c2′lnt+c2′=t6+34⟹c2′=t6+34
回代求 c1′:
c1′=−(t6+34)lnt=−t6lnt−34lnt
Step 4: 积分
- c2=6lnt+34t
- c1=−3ln2t−34(tlnt−t) (分部积分)
Step 5: 组合通解
将 c1,c2 代入 x=c1t+c2tlnt 整理即可。
题目意思:
现在的基本解组 et,e−t 在 t=0 时的数值乱七八糟的。
我们想要两个新的“完美”解 ϕ1(t) 和 ϕ2(t),它们满足:
- ϕ1(0)=1,ϕ1′(0)=0 (这叫“位置解”,只管初始位置)
- ϕ2(0)=0,ϕ2′(0)=1 (这叫“速度解”,只管初始速度)
解法:
第一步:构造 ϕ1
设 ϕ1(t)=c1et+c2e−t。
代入初值条件:
[
{c1+c2=1(函数值)c1−c2=0(导数值)
]
解得 c1=21,c2=21。
所以 ϕ1(t)=2et+e−t=cosht (双曲余弦)。
第二步:构造 ϕ2
设 ϕ2(t)=k1et+k2e−t。
代入初值条件:
[
{k1+k2=0k1−k2=1
]
解得 k1=21,k2=−21。
所以 ϕ2(t)=2et−e−t=sinht (双曲正弦)。
第三步:利用标准解组写任意解
有了这两个标准解,任意初值 x(0)=x0,x′(0)=x0′ 的解直接写出来,不用计算:
x(t)=x0⋅ϕ1(t)+x0′⋅ϕ2(t)x(t)=x0cosht+x0′sinht
这是一个理论证明题,核心在于行列式怎么求导。
证明思路:
行列式求导法则:等于“每次只对一行求导,其他行不变”的行列式之和。
对于 n 阶朗斯基行列式 W(t):
- 第一行是 xi,第二行是 xi′,……,最后一行是 xi(n−1)。
- 对 W(t) 求导,会得到 n 个行列式的和。
- 关键点:前 n−1 个行列式中,都会出现两行完全一样的情况(比如对第一行求导变成 xi′,这就和原来的第二行一样了)。两行一样的行列式等于 0。
- 唯一剩下的:只有对最后一行求导的那个行列式不为0。
[
W′=x1⋮x1(n−2)x1(n)⋯⋯⋯xn⋮xn(n−2)xn(n)
]
- 代入方程:因为 xi 是方程 x(n)+a1x(n−1)+⋯+anx=0 的解,所以:
xi(n)=−a1xi(n−1)−⋯−anxi
- 把这个代入到 W′ 的最后一行。利用行列式拆分性质,只有包含 −a1xi(n−1) 的那一项是非零的(因为其他的项如 −anxi 会让最后一行和第一行成比例,行列式为0)。
- 提取公因子 −a1,剩下的正是原行列式 W(t)。
- 结论:W′=−a1(t)W。这是一个简单的一阶方程,解它就得到公式。
直觉理解:
想象三维空间(n=3)里的一个平面。
- 齐次方程的解通过原点,构成一个平面(3维线性空间)。
- 非齐次方程的解,就是这个平面平移了一下(不过原点)。
- 在这个平移后的平面上找点,你最多能找几个点,使得它们不“共面/共线”?
证明逻辑:
- 存在性:取齐次方程的 n 个线性无关解 x1,…,xn 和非齐次方程的一个特解 xˉ。
构造 y0=xˉ, y1=xˉ+x1, …, yn=xˉ+xn。
这 n+1 个解是线性无关的。(可以用定义证)。
- 最多 n+1 个:假设有 n+2 个解 y1,y2,…,yn+2。
考虑它们的差:z1=y2−y1, z2=y3−y1, ….
根据非齐次解的性质2,这些差 zi 都是齐次方程的解。
齐次方程的解空间维度是 n。
如果你有 n+1 个差向量(由 n+2 个解产生),它们在 n 维空间里必然线性相关。
既然差向量线性相关,那么原解 yi 也就线性相关了。
所以,不可能有 n+2 个线性无关解。
希望这套详细的讲解能帮你打通任督二脉!如果哪个积分算不出来,或者哪个步骤卡住了,继续问!